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Abstract 

Time-dependent correlation functions of (unstable) particles undergoing biased or un- 
biased diffusion, coagulation and annihilation are calculated. This is achieved by similar- 
ity transformations between different stochastic models and between stochastic and soluble 
non- stochastic models. The results agree with experiments on one-dimensional annihilation- 
coagulation processes. 



PACS: 02.50.Ga, 05.40. +j, 82.20.Mj 



The physics of interacting particles out of thermodynamic equihbrium has been a subject of 
much recent interest. While in larger spatial dimensions, conventional rate equation approaches 
are sufficient, systems constrained to be effectively one-dimensional display novel and interesting 
fluctuation effects. For example, for particles A diffusing on a lattice and undergoing a binary 
reaction process A + A ^ products one expects for large times t an algebraic fall-off of the mean 
particle concentration 

c{t) ~ ry (1) 

In ID systems, one finds ?/ = | as opposed to y = 1 obtained from a (mean-field) rate equation. 
The exponent y has also been measured experimentally in effectively one-dimensional systems. For 
annihilation-coagulation reactions A + A ^ products one finds y = 0.52 — 0.59 and y = 0.47(3) 
[0, and for the (pure) coagulation reaction A + A — A, ?/ ~ 0.48 0. Although these reacting 
systems might appear to be quite different, in this work we show that these and more general 
systems can be treated in a simple and unified way. In particular, a simple explanation for the 
same value of y in all annihilation-coagulation problems is obtained. 

For the theoretical description of these reaction-diffusion systems, a useful approach g] consists 
of rewriting the master equation which describes the time evolution of the probability distribution 
function P{{f3};t) Schrodinger equation 

dtPi{/3};t) = -HPm;t) (2) 

in which the quantum Hamiltonian H is defined in terms of the transition rates w{(3 ^ '-/) between 
two configurations {j3} and {7} by 

(7| H = -w{[3 ^ 7) , H\I3) = Y: w{(3 ^ 7) • (3) 

H describes a stochastic process since the sums of all elements in each column add up to zero. 
This conservation of probability is equivalent to the relation 

{s\H = , {s\=Y.m\ (4) 

for the vector {s\. Then the well-known theorems about the solutions of the master equation [^] 
can be translated into the Hamiltonian formulation at hand. In particular, the real parts of the 
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eigenvalues of H are non-negative. Furthermore, starting from an initial probability distribution 
\Po) = J2i3 P{{P}j t = 0) where each configuration {/?} occurs with probability P({/?}; t = 0), 
the solution to the master equation is then formally given by the time-dependent probability 
distribution (state vector) \P) = Y.f3 P{{P}jt) = exp{—Ht) \Po)- Time-dependent averages 
of an observable^ F are given by the matrix element 



<F>(t) = (s|Fe-^*|Po) (5) 

The interest in this setup of the problem in ID comes from the integrability of the quantum 
Hamiltonian H for large classes of reaction-diffusion processes P, 0) §• Exactly known results 
such as the knowledge of the spectrum of H (obtained e.g. from the Bethe ansatz) have so far 
led to a number of exact and explicit results for time-dependent averages and correlations. Here 
we show how these and other existing results (see e.g. ^, |T0| for annihilation and coagulation 
processes) can be extended to considerably more general quantum Hamiltonians: 

1. We investigate relations between stochastic systems whose quantum Hamiltonians H and H 
are related through a change of basis of the one-particle states, see |^, ^, |Tl], |12 



H = BHB-^ , B = ^Bi (6) 



i=l 



where Bi is acting only on the site i. 

2. We consider a given quantum (and in general non-stochastic) Hamiltonian with known 
properties and we look for stochastic processes which can be obtained from this Hamiltonian 
by a similarity transformation of the form (P). The time- dependent behaviour of these new 
stochastic systems can then be elucidated in terms of the original Hamiltonian. 

Under the similarity transformation B from (P) averages transform as follows 

<F>{t) = {s\ Fe-^' \Po) = {s\ Fe-^' \Po) (7) 
with the transformed observable F = FB~^ and transformed initial distribution Pq\ = B \Po)- 



is a suitably chosen projection operator (see below). 
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We now give the general form of the single-site matrix B for transformations between two 
stochastic systems. Certainly, one-site state vectors of the system S described by the Hamiltonian 
H must have the form |p) = ^ ^ ^ ^, with < p < 1. Also, for the transformed state B |p), 
probabilities must sum up to one, for all values of p, thus 

bu{l -p)+ &21P + buil -p)+ &22P = 1 (8) 
where the bij are the elements of B. Comparing coefficients, we get for B the form 

Since obviously {s\ B = {s\, it follows that {s\H = for the full system defined on L sites. 

The stochastic systems S = {H, p) under consideration are described by the quantum Hamil- 
tonian H and a set of parameters p specifying the initial conditions (see below). The relations 
between two stochastic systems 5* and 5* are caught by the following two definitions. 

1. The transformation S —>■ S between two stochastic systems S and S is called a similarity 
transformation, if there exists a non-singular transformation B of the form between the 
quantum Hamiltonians H and H such that all reaction-diffusion rates are positive in both 
systems. S and S are then called similar. 

2. The transformation S ^ S between two systems S and S is called a stochastic similarity 
transformation (SST), if S and S are similar and furthermore if for all probability distribu- 
tions \P) oi S also P) = B \P) is a probability distribution of S. 



Through similarities and SST a given system may be mapped into a simpler form. Examples will 
be given below. To illustrate the second definition, consider an uncorrelated initial state of the 
form 

iPo) ^ g ( ^ ) (10, 

with < Pi < 1 for all sites i = 1, . . . , L. If 5* is obtained from 5 by a SST, initial states of the form 
(|10D are mapped onto transformed initial states Pq) of the same form and with < p^ < 1 for all 
sites. We stress that the notion of a SST between two stochastic systems 5* and 5* is considerably 
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diffusion to tlie left 
diffusion to the right 
pair annihilation 
coagulation to the right 
coagulation to the left 
death 

decoagulation to the left 
decoagulation to the right 
birth 
creation 
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Table 1: Two-sites reaction-diffusion processes and their rates. 



more restrictive than mere similarity, which does not require that also the transformed state vector 
is a probability distribution. We also remark that because of the locality of the change of basis 
the results obtained here are valid in an arbitrary number of space dimensions, although we shall 
present the argument explicitly only for (i = 1. 

We now state precisely the models we shall study below. Consider a one- dimensional lattice, 
with L sites and periodic boundary conditions. Each lattice site can either be empty (denoted by 
0) or occupied by a single particle (denoted by A). Particles can hop to an empty nearest neighbor 
site. A single particle or a pair of particles on neighboring sites can undergo a chemical reaction. 
The reactions we are going to consider are specified with their rates in table |1]. 

Furthermore, we shall distinguish between unbiased and biased reactions. For unbiased re- 
actions, the rates with indices L and R are all equal and we shall then drop the index (e.g. 
Di^ = Dr = D etc.). For biased reactions, we define an anisotropy parameter rj from 

DL = D{l + r]) , 7i = 7(l + r/) , /5^ = /5(l+r/) 

DR = D{l-r]) , 7^ = 7(1-,^) , /5« = /5(l-r/) (11) 

For ?7 = we recover the unbiased case. 

We begin with unbiased systems, that is r/ = 0. We consider the following system, with 
diffusion, annihilation, coagulation and death reactions present (see table |I|). In what follows, we 
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take units of time such that D = 1. Then the Hamiltonian can be written in terms of two-site 
contributions H = Y^i^i^i+i, where 



/ 


-5 


-S -2a \ 





1 + 5 


-1 -7 





-1 


1 + 6 -7 


V 





2(a + 7) / 



(12) 



We define the A;— point correlation function of the particle number operator n^^., i = 1, 2, . . . , /c as 

Ck{{x},t;a,-f,6; p) =< n^^ ■ ■ -n^.^ >h (t) (13) 

where we explicitly indicate the dependence on the rates as well as on the initial conditions. The 
operator rix is a projector with eigenvalue 1 if site x is occupied and eigenvalue if it is vacant. 
Although H is non-hermitian, it is known that there is a decomposition H = Hxxz + Hi into an 
hermitian Hamiltonian Hxxz (which is the Hamiltonian of the anisotropic Heisenberg quantum 
spin model) and a non-hermitian part Hi such that the eigenvalues of H are exactly the eigenvalues 
of Hxxz @- That is so because the chemical reactions permitted here only destroy and never 
create particles. In one dimension, an interesting special case is given by the free-fermion condition 

a 7 = 1 + (5. (14) 

In that case the hermitian part Hxxz can be diagonalized exactly in terms of free fermions. If 
either just annihilation or coagulation are present, it is known that a closed system of equations 
of motion can be found 0. (|1^ means that diffusion and death together occur at the same rate 
as annihilation and coagulation together. If 5 < 7, we can rewrite the problem as an annihilation- 
coagulation problem of an unstable particle, where the effective coagulation rate is modified into 
7efr =1 — 6, and 1/(25) is the life time of the unstable particle.0 If we use the diffusion process to 
determine the time scale, we can say that if two particles attempt at the same time to be on the 
same site, they undergo a chemical reaction with probability one. The ratio jes/ct then gives the 
branching ratio between the two processes. 

At first sight, the condition ( |I4D appears to be rather artificial. However, it is apparently 

realised to good approximation in one of the experimental realizations of the model considered 
■^The special case 7 = 5, i.e. pair annihilation of unstable particles is discussed in the appendix. 
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so far 0]. The carrier substance is N{CH-i)4^MnCl^. The particles are excitons of the Mn"^^ ion 
and move along the widely separated MnCl^ chains. A single exciton has a decay time of about 
0.7ms. The on-chain hopping rate is 10^^ — lO^^s^^. If two excitons arrive on the same Mn^^ ion, 
they undergo a coagulation reaction with a reaction time ~ 100/s Since the reaction time is 
much smaller than the diffusion time, we can conclude that the reaction probability is very close 
to one. This gives back (]T^), with a = 5 = for this example. 



After these preparatory remarks we return to the general case. Starting from the system S as 
defined through its Hamiltonian ([1^), we get the following simplified systems S through a SST. 



I. (5 > 2q; + 7. Through a SST we get the system Si with 

^ ^ ~ ^5 — Icy. — 7 

D = D , q; = , 7 = q; + 7 , 5 = 5 , p = -j p (15) 

— a — 7 

and we find 

(5 — ol — 7^^^^ 5 — lot — 7 
T — 7, C'fc({a;},t;0,a + 7,5;- p) (16) 

— za — 7/ — a — 7 

II. (5 = 2a + 7. In this case, the transformation becomes singular. However, the equations 
of motion for the particle number correlators C({a;}) decouple from each other ^j. For 
example, the particle density at time t is 

00 L 

Ci(a;,t)=e-2(i+^)* ^ ^ Ci(y, 0)W^^(2t) (17) 

m=— 00 y=l 

where is a modified Bessel function. 
III. (5 < 2a + 7. We find the transformed system Sni 

D = D , a = a + 7 , 7 = , 5 = 5 , p = ^^±2^ p (18) 

2a + 27 — d 

and we find 

/2a + 27-5\^ 2a + 7 -5 

Cki{x}, t;a,j,5;p) = ( + ^ _ ,5 ) '^fc({a^}^ t; a + 7, 0, ^'^ 2a + 2^-5 ^'^^^ 
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IV. (5 < a + 7- In this interval, there is a second SST onto the system Siy with 

D = D , 5 = a + 7-5 , 7 = 5 , 6 = 5 , p= ^ p (20) 

2a + 27 — 26 

and we find 

Cfc({x},t;a,7,(5;p) = ( — ^) C^do:}, t; a + 7 - 5, (5, (5. ^ — ^p) (21) 

\ za + 7 — y 2(a; + 7 — 0) 

Although both systems 5*///, Siv are found from 5* through a SST, the latter is more useful for 
practical calculations. Using the results derived in the appendix, we can isolate the dependence 
on 6 completely. For translationally invariant initial distributions we have for the large-time 
behaviour in the free fermion case (see the appendix for more general cases) 

Ci(x,t) - poe-''* (22) 
C2ix,x + r,t) ~ r3/2g-45t . ^2^^ 

and we explicitly see that for 6^0 the initial particle density does enter into the large-time 
behaviour. 

This is different from the result found when the death reaction is absent {6 = 0). In that 
case only the SST onto Sjjj = Sjv remains. (The corresponding similarity transformations have 



been derived before [jTT], [T^.) All correlations then depend non-trivially only on the effective 
annihilation rate Ofeg = + 7. The correlation function C on the r.h.s. of ( PT| ) then is the known 
density correlation function for diffusion annihilation, see ^ 0]. For example, with the initial 
state (p!OD with Pi = p for all sites i we have for the mean particle concentration c{t) ~ / dxCi{x, t) 
in the free fermion case for the process A + A ^ 



c{t) = pe 



-ADt 



MDt)-^'^ (1 + 0{t-^)) (23) 



/o(4Dt) + 2(1 - p) ^(1 - 2pf-^h{ADt) 

k=l 

In particular, we always get back ?/ = ^ in (|l]), in agreement with experiment [IH @, 0] • Furthermore, 
the data of Kroon et al. show that the long-time behaviour of c{t) is independent of the initial 
particle density p, in agreement with (^). 

So far, the transformations considered have mapped S back onto itself, up to changed values of 
its parameters. But it is sometimes possible to reduce more complex systems to the ones discussed 



here. For example, the system S with the parameters 



D = l , (5 = 2a + 27 (24) 

and 7 7^ gives through a SST the system 5* with (see table |l]) 

D = \{3 + 2a + 2-f) , 5 = , 7 = ^(0 + 7) , 5 = 27 , V = 4^ (25) 

Since the system contains particle creation as well as particle destruction terms, it no longer 
has a trivial (i.e. empty) steady state. This steady state can be found easily, since for a single-site 

state 

Since the steady state of 5* is just {8> ^ g ^ , we find that the steady state of S has a mean particle 
density p = 2/3. The approach towards this steady state is exponential. 

Another example is found when a + < 6 < 2a + 6. Then 5* is similar to 5* with 

D = ^(5-a-7) + l , fi = l(5-a-7) , ^ = s = a = ^{2a + 2^ - 6) , ^=^6 (27) 



and the one-site state changes into 



and we get a steady state particle density of p = 1/3. The approach towards this steady state is 
exponential. The transformation 5* ^ S* is a SST if 5 > |a + 7. Conversely, 5* — > S* is a SST if 
5 < |a + 7. Other examples with 6 = are given in |I2| . 

We now turn our attention to some systems with biased react ion- diffusion processes. We take 



diffusion, coagulation and annihilation into account. The rates are given in (|TT]). Using the 
unbiased case 17 = as a guide, we seek a SST S ^ S such that 7^ = 0. In fact, using 61 = 1 as 
before and choosing 62 in order to get 7/^ = 0, we find 

A(a + 7 — 1)7^7 _ ~ , , 

7L = ^ ! , V = 7^ , D = D = 1 29 

2a + 'j + 77(7 - 2) 
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For = 0, we recover the previous result. However, if we use the free fermion condition a + 7 = 1, 
then 7l = and 5 = 1. We then have 



Ck{{x},t; a, 'J 



a-rj-p) 



l + a 



C,({x},t;l,0;r/;^(l + a)p) 



(30) 



Considering the mean particle density only, this relation was also observed in for the special 
case of mapping the pure biased coagulation problem (a = 0) onto the pure biased annihilation 
problem (7 = 0). Generally one finds that for a homogeneous initial condition (10) with pi = p = 
const the correlation function Ck{{x}, t; 1, 0; 77; (1 + a)p/2) is independent of the bias r] [|10]. For 
an inhomogeneous initial state with = 1 and Py = 1/2 ioi y ^ Xq one finds for the average 
density in an infinite system for large times 



Ci{x,t) =< > (t) 



1 



+ 



-(x-xo-r/t)'^ /2t 



(31) 



/47rt 7rt2 

Finally, we illustrate the transformation (|^) between a stochastic and a non-stochastic system. 
As an example, consider the Hamiltonian H = J2i 



i,i+l 



f A -2a 

A-1 -D 

-D A-1 

V A-2 J 



(32) 



H can be solved in terms of free fermions. We want the transformed Hamiltonian H to describe 

( }) }) \ 

a stochastic system, that is we require that {s\H = 0, see (H). Writing B = \ , , , the 
solution to this condition is 



J21 



-h 
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A 



(33) 



We now take 612 = 622- Let F = (611/622)^ > 0. Then the positivity of the reaction rates in H 
requires that T = D/a. The Hamiltonian then reads 



.i+l 



/ 2-2D -1 - D -1 - D \ 

-1 + D 2 + 2D D-1 

-1 + D 2 + 2D D-1 

V -1- D -1- D 2-2D J 



(34) 



The off-diagonal elements of H are non-positive provided < D < 1 and < a. Under these 
conditions B is non-singular and H is the quantum Hamiltonian of a stochastic system. We point 



out that this Hamiltonian is identical to the quantum Hamiltonian for the ID Glauber-Ising model 
[0 



Hgi = E (1 - (l - f (o^- c^f+i + ^li^t)) (35) 

at temperature T given hj D = tanh (2 J/Zc^T). In this way we obtain a new relation between 
non-zero temperature Glauber dynamics and the XXZ chain in the free fermion case. On the 
other hand, for non-vanishing T the Glauber-Ising model can be transformed into an XY free 



fermion chain [|T5|. We shall return to consequences of this observation and the reformulation of 
more general stochastic processes in terms of soluble free fermion systems elsewhere. 

In summary, we have shown how to relate different stochastic systems using similarity transfor- 
mations. In several examples, this technique proves to be useful to extend the scope of integrable 
systems. The results are in agreement with the existing experiments and include some previ- 
ous observations of relations between different systems as special cases. Going beyond similarity 
transformations between stochastic systems, we have found a simple example how to reformu- 
late a stochastic system in terms of a non-stochastic soluble free-fermion model in a novel way. 
The techniques developed here can be used for a systematic study and classification of stochastic 
interacting particle systems in terms of integrable quantum chains. 
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Appendix 



Here we study the annihilation process for unstable particles with an average life time 7 = 5. The 
equations of motion for the average density and density correlations read 
d 

— < > = < n^.+i > + < n^-i > -2(1 + d) < > 

-2a (< n^n^+i > + < n^-iu^ >) (A.l) 
^ < n^Uy > = < n^+iUy > + < n^-iHy > + < n^riy+i > + < n^Uy-i > -4(1 + 5) < n^Uy > 
-2a (< n^nr^+iHy > + < rir^-irir^ny > + < n^UyUy+i > + < n^ny-iUy >) 

(if|?/-x|>2) (A.2) 
< n^Ux+i > = < n^-in^+i > + < > -2(1 + a + 25) < n^n^.+i > 



dt 



-2a (< n^-in^n^+i > + < >) (A. 3) 



and similar expressions for higher order correlators P|. For a A;-point correlator there is always 
a coupling to {k + l)-point correlation functions proportional to the annihilation rate a. First 
we show that the amplitude of a fc-point density correlation function Cfc({x}; t) =< n^^ . . . n^^ > 
decays for large times t at least with a factor proportional to exp {—2kSt). 

To see this, we recall that the spectrum of the quantum Hamiltonian H is exactly the same 
as of the XXZ quantum chain ||^. For the calculation of the eigenvalues, it is thus sufficient 
to consider the sectors with fixed number of particles k separately. In the fc-particle sector, the 
eigenvalues of H are 

N 

= fc2(5 + ^2(l-cosgi) > A;25 (A.4) 

i=l 

and the values of the Qi are determined from the Bethe ansatz equations. On the other hand, the 
A;-point correlators can only take non-zero values when defined on states which contain at least 
k particles. Furthermore, Ck depends through the equations of motion directly only on Ck and 
Cfc+i and in particular it is independent of Cq, Ci, . . . , Ck-i- Thus, for large times t we must have 

Cki{x};t) ~ e"^* , with A > k25 (A.5) 



because of the inequality ( |A.4| ) for the eigenvalues of H in the A;-particle sectoiQ. 



^Writing Ck{{x};t) = e ^''^'^Bk{{x};t), one can further show with the Bethe ansatz that genericaUy Bi{t) 
0(1) and Bk>2{t) — > as i ^ oo. 
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Now we define the explicitly time-dependent quantity n^it) = exp {2St) n^^. According to 
the considerations of the previous paragraph A;-point correlation functions of this quantity are 
bounded from above by their initial values < n^.^ . . . h^^ > (t = 0) =< . . . n^^ > (t = 0) < 1. 
Rewriting Eqs. ( [A.ID - ( |A.3| ) and the corresponding equations for /c-point correlators in terms 
of averages for n^- reproduces equations of the same form, but with effective time-dependent 
couplings a(t) = aexp {—26t) to {k + l)-point correlators and with 6 = 0: 

d ^ _ . ^ , 

-— < > = < n^^i > + < n^-i > —2 < > 



dt 



-2ae (< n^n^+i > + < n^-in^ >) (A.6) 
— < n,^ny > = < n^+iUy > + < n^-iUy > + < n^Uy+i > + < n^Uy^i > -4 < n^Uy > 

-2ae"^'^* (< n^na;+iriy > + < h^-ih^hy > + < h^fiyhy+i > + < h^hy-ihy >) 
(if \y-x\> 2) (A.7) 
— < n^n^+i > = < n^_in^+i > + < > -2(1 + a) < n^n^+i > 

-2ae"^'^* (< fix-in^n^r+i > + < r2^n^.+ifi^+2 >) (A.8) 

Since < . . . hx^. > (t) < 1 for all times t and since for long times ae"^''* can be neglected 
these equations effectively decouple and reduce to closed linear differential-difference equations 
for /c-point correlators which can be solved with the Bethe ansatz [Q. 

The equation for the one-point function reduces to a lattice diffusion equation which is solved 
by modified Bessel functions. Assuming a translationally invariant initial state with 
< rixit = 0) >=< n^it = 0) >= po one finds for the average particle density at time t 

< n,{t) >^ poe-^'' (A.9) 

This quantity depends on the initial density which is in contrast to the diffusion limited anni- 
hilation of stable particles where the density decays algebraically (for long times) and with an 
amplitude independent of the initial density. 

Defining C{r,t) =< hxfix+rit) > and using ( |A.6| ) - ( |A.8D one finds for the two-point function 



with translationally invariant initial conditions the equations 

d_ 
di 



^ Cir,t) = 2{C{r + l,t) + C{r-l,t)-2C{r,t)) (r > 2) 
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jC{l,t) = 2(C{2,t)-{l + a)C{r,t)) (A.IO) 

This gives for the two-point density correlation function C{r,t) —< nxUx+rit) > in an infinite 

system 

oo 

C(r, t) ^ e-<^+'5)* Yl K^r-ym + bylr+y-im] (A-H) 
y=l 

where In is the modified Bessel function, ay are constants defined by the initial distribution and 
by = fiay — (1 — /U^) ^21=1 IJ^^~^~^cik with n = 1 — a. At first sight, we should expect for large times 
that C{r,t) ~ exp(— 45^). For the free fermion case a — 1, however, a different result is 
found. Since by — —Uy-i we get 

oo 

C{r, t) « e-^(^+^)* J2 [^-.(4^) - ^.+.(4^)] (A.12) 

y=i 

For <^ t this correlator decays as C ~ exp {—A5t) whereas for r"^ t one has C ~ 
exp {—ASt). The same effect is also seen for a — 2. 
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